A stacking-fault based microscopic model for platelets in diamond 
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We propose a new microscopic model for the {001} planar defects in diamond commonly called 
platelets. This model is based on the formation of a metastable stacking fault, which can occur 
because of the ability of carbon to stabilize in different bonding configurations. In our model the 
core of the planar defect is basically a double layer of three-fold coordinated sp 2 carbon atoms 
embedded in the common sp 3 diamond structure. The properties of the model were determined 
using ab initio total energy calculations. All significant experimental signatures attributed to the 
platelets, namely, the lattice displacement along the [001] direction, the asymmetry between the 
[110] and the [110] directions, the infrared absorption peak B' , and broad luminescence lines that 
indicate the introduction of levels in the band gap, are naturally accounted for in our model. The 
model is also very appealing from the point of view of kinetics, since naturally occurring shearing 
processes will lead to the formation of the metastable fault. 

PACS numbers: 61.72.Nn, 61.72.Bb, 63.20.Pw, 71.15.Nc, 71.55.Cn 



Perhaps one of the oldest unanswered questions in di- 
amond physics concerns the nature and atomic struc- 
ture of the extended defects known as platelets. These 
{001} planar defects were first discovered more than sixty 
years ago by X-ray diffraction experiments, which ob- 
served anomalous peaks corresponding to {00ft.} reflec- 
tions. These were immediately associated with lattice 
defects, since they are forbidden by symmetry in a per- 
fect diamond lattice. From then on, a plethora of experi- 
mental data on this planar defect has been gathered, but 
a complete understanding of its origin and microscopic 
structure is still lacking. The current knowledge on this 
defect can be summarized as follows: i) platelets have 
been detected only in type la diamonds |2| ; ii) transmis- 
sion electron microscopy (TEM) experiments have shown 
that there is an asymmetry between the [110] and [110] 
directions 3]; hi) TEM experiments have also determined 
that platelets displace the crystalline lattice by approxi- 
mately 0.4ao [4|, along the [001] direction, where a is 
the lattice parameter of diamond; iv) electron energy 
loss spectroscopy (EELS) experiments indicate that the 
nitrogen content inplatelets can vary from 6% to 61% 
of a monolayer 0, @, ; v ) & n infrared absorption line 
around 1370 cm -1 is always present in samples contain- 
ing platelets 0; vi) platelets have been associated with 
broad luminescence bands, one centered at 1.25 eV 0, 
which reduces the efficiency of optical windows made of 
natural diamond, and another at 2.14 eV [10J, and possi- 
bly with high energy absorption and luminescence bands 
around 4.6 and 4.4 eV, respectively |ll| . 

In the microscopic model first proposed for platelets 
in diamond, Frank 01 considered that the defect should 
be formed by Si impurities replacing part of the carbon 



atoms in a (001) plane. However, it was later shown that 
Si is an uncommon impurity in diamond. Subsequently, 
it was demonstrated that diamonds type la contain N 
impurities in the required concentration [13j. Follow- 
ing that, Lang |l4| proposed that the defect would be 
constituted by a double layer of N atoms (one layer of 
substitutional N atoms and one layer of interstitial N 
atoms). However, as mentioned above, the N content 
in the platelets can vary substantially 0, E3, 0i ruling 
out N impurities as the main constituent of the defect. 
Since foreign impurities could not account for the ex- 
istence of platelets, Evans 15] proposed that platelets 
could be formed by interstitial carbon atoms. Based on 
Evans' idea, Humble suggested that the double layer 
of nitrogen atoms in Lang's model should be replaced by 
a double layer of carbon atoms, thus establishing a so- 
called interstitial-aggregate model that has since become 
the most accepted microscopic model for platelets in dia- 
mond. More recently, Baker proposed another model 
for platelets based on the aggregation of the Rl centers. 

Very recently, several interstitial models, based on 
Humble's proposal, have been extensively studied by 
Goss and co-workers |l7j . using ab initio calculations. 
Their results indicate that some of the experimental sig- 
natures of the platelets can be explained by the intersti- 
tial model. However, the aggregate of interstitials does 
not introduce states in the electronic band gap, and hence 
it alone cannot account for the optical activity experi- 
mentally observed. Goss et al. propose that the observed 
levels in the gap are due to vacancies and nitrogen im- 
purities segregated at the platelets. Therefore, one could 
expect these properties to be somehow sample depen- 
dent. In fact, the main difficulty with the interstitial 
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aggregate model concerns the energetics and the atom- 
istic mechanism of its formation. The formation energy 
of interstitial defects in the bulk is quite high (~ 12 eV). 
Goss et al. have proposed that forming a Frenkel 

pair (a complex of an interstitial and a vacancy) at the 
platelet, would result in a substantial reduction in the 
interstitial formation energy. 

Here we propose a new model for the microscopic struc- 
ture of platelets in diamond based on a entirely different 
mechanism. The model does not require the creation of 
point defects, instead, the defect is created by a shear- 
ing process of the crystal. It sprung from the study of 
the 7-surface 0] or the generalized stacking fault en- 
ergy landscape for the (001) slip plane in diamond. In 
particular, we studied the behavior of the 7-surface along 
the [110] and [110] directions using ab initio calculations. 
The model, aside from accounting for all the known ex- 
perimental signatures of platelets in diamond, requires a 
much smaller activation energy for the formation of the 
planar defect than the interstitial aggregate model. 

The results we present here were obtained through 
ab initio calculations, within the framework of the den- 
sity functional theory (DFT) and the generalized gradi- 
ent approximation (GGA), using the SIESTA code [l9j . 
A double-zeta localized basis set with polarization or- 
bitals was used. The interaction between the valence 
electrons and the ionic cores was modeled using ab ini- 
tio norm-conserving pseudopotentials. A 216-atom su- 
percell was used in most of the calculations, except for 
the vibrational density of states calculations, where a 64- 
atom supercell was employed. The T point was used in 
the Brillouin zone sampling for the 216-atom supercell, 
whereas for the calculations using the 64-atom supercell 
8 k-points were used. To check the accuracy, we calcu- 
lated the intrinsic stacking fault energy in the (111) slip 
plane, obtaining 0.279 J/m 2 , in excellent agreement with 
the experimental value 0.285 ± 0.040 J/m 2 |2c| . 

The 7-surface for a given crystalline plane is obtained 
by cutting the infinite crystal in half along the plane and 
shearing the upper part with respect to the lower part by 
a displacement (slip) vector, which belongs to the plane. 
Our results for the cut of the 7-surface along the [110] 
direction are shown in Fig. 1. One of the curves displays 
the unrclaxcd 7-surface, which exhibits a maximum cor- 
responding to a displacement of (arj/4)[110]. For this dis- 
placement, the atoms in both sides of the cut approach 
each other very closely and are aligned along the [001] 
direction. Fig. 1 shows a quite large release of energy 
that results from allowing the atoms to relax along the 
[001] direction only. In fact, we found that the slip at 
(oo/4)[110] corresponds to a local minimum of the 7- 
surface along this direction, in contrast with the case of 
Si and Ge, where this configuration corresponds to the 
so-called unstable stacking fault. On the other hand, 
shearing along the [110] direction does not bring the sys- 
tem to any intermediate stable configuration, and our 
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FIG. 1: Profile of the 7-surface of the (011) plane along the 
[110]. The inset shows the details of the energy minimum for 
the displacement (ao/4)(110) after atomic relaxation. 



calculations confirm that the 7-surface does not exhibit 
any local minimum, a behavior in this case similar to 
what is found in silicon and germanium |2*lj . Since at 
first we allowed relaxation of the (arj/4)[110] configura- 
tions along the [001] direction only, to verify that this 
is a local minimum of the 7-surface we allowed all the 
atoms to move freely. We found that, indeed, the atomic 
structure is stable. That such configuration is stable in 
diamond and not in silicon and germanium, stems from 
the stability of the sp 2 bonding in diamond. 

Fig. 2a depicts the relaxed atomic structure corre- 
sponding to the metastable stacking fault at (ao/4)[ll0]. 
The blue colored carbon atoms are three-fold coordi- 
nated, whereas the gold colored ones are the usual four- 
fold coordinated carbon atoms in diamond. The core of 
the defect is formed by a double layer of three-fold co- 
ordinated carbon atoms. In Fig. 2(b), we display the 
structure of one unit of the core of the defect (with a few 
neighboring atoms), when viewed along the [110] direc- 
tion. This unit is formed by the repetition of a hexagonal 
ring consisting of four three-fold coordinated atoms and 
two four- fold coordinated atoms. The length of the bonds 
between three- fold coordinated atoms is 1.35 A, which is 
even shorter than the carbon bonds in graphite (1.42 A), 
while the bonds between three-fold coordinated and four- 
fold coordinated atoms are about 1.52 A long, very close 
to the value of 1.57 A we obtained for the bond length 
in bulk diamond. Note that, in Fig. 2(b), the length 
of the bonds between the four-fold coordinated atoms 
just outside the defect core is already the bond length 
in bulk diamond. Therefore, this structure can be seen 
as two semi- infinite pieces of diamond "glued" together 
by chemical bonds that are 1.35 A long, which is equal 
to 0.38cto, basically the lattice displacement observed in 
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FIG. 3: Vibrational density of states in bulk diamond and in 
diamond containing a planar defect 




FIG. 2: (a) Atomic structure of the planar defect. The 
blue colored atoms are three-fold coordinated carbon atoms 
and the gold colored atoms are four-fold coordinated carbon 
atoms, (b) Detail of the core of the planar defect indicating 
bond lengths and angles. 



platelets. Moreover, it is also easy to see that the atomic 
structure of the core along the [110] direction is different 
from that along [110] direction, reflecting the asymme- 
try of the 7-surface between the two directions. This 
extended defect has been discussed before [2^, however, 
within a completely different context, related to grain 
boundaries in artificially grown diamonds. 

In Fig. 3, we show our ab initio results for the vi- 
brational density of states (VDOS) for the planar defect 
and for bulk diamond. First, we remark that our calcu- 
lations can reproduce quite well the bulk VDOS. More- 
over, our results show clearly the appearance of three 



peaks above the highest allowed frequency in the bulk. 
Two of these peaks lie just above the bulk band edge, 
at 1373 cm -1 and 1441 cm -1 , values which are in very 
good correspondence with the frequencies (1372 cm -1 
and 1426 cm -1 ) of the experimental bands usually asso- 
ciated with platelets[8(. The first peak at 1373 cm -1 is 
more intense and can be associated with the B' band. By 
computing the vibration frequency of the stretch mode 
of the compressed bonds in the core of the defect, we 
determined that the third peak in Fig. 3, at 1854 cm -1 , 
is a localized mode associated with these compressed sp 2 
bonds in the core [24| . As far as we know this frequency is 
out of the range usually investigated in infrared absorp- 
tion studies of platelets, 900 cm -1 to 1650 cm -1 [23j. 
This is quite interesting, because this third peak opens 
the possibility that the validity of our model could be 
experimentally investigated. 

We now discuss the electronic structure of the 
metastable fault. Fig. 4 depicts both the electronic den- 
sity of states of a bulk diamond cell and of a supercell 
containing the planar defect. In the figure, the zero in the 
energy scale corresponds to the highest occupied level. 
The first important point to note is that the planar de- 
fect gives rise to deep levels in the band gap. The inset in 
Fig. 4 shows the levels in the diamond band gap without 
broadening. The lowest empty levels lie about 1.0 eV. 
Also, in the inset, one can see empty levels about 2.1 eV. 
Above that, there are empty levels that coincide with 
the bottom of the conduction band of the bulk at about 
3.1 eV (not shown in the inset). Our results agree with 
the luminescence bands at 1.25 eV and 2.14 eV, usually 
associated with the platelets. It is important to empha- 
size that these levels appear in the gap as consequence 
of the planar defect, without the intervention of either 
native defects or impurities. 

We have also investigated the possibility of the oc- 
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FIG. 4: Electronic density of states of bulk diamond and of 
diamond containing a planar defect. The inset shows the 
levels without broadening. 



currence of a similar intermediate metastable stacking 
fault in the (111) slip plane. For this plane, a very sim- 
ilar situation, with the atoms adjacent to the cut get- 
ting near each other, also happens when the system is 
sheared along the [121] direction by a displacement of 
(\/6/12)ao. In the case of Si and Ge, previous calcula- 
tions |2l| indicate that the atomic configuration obtained 
after relaxation is still a distinct local maximum for the 
7-surface of the glide set. In the case of diamond, the 
7-surface resembles a plateau, which turned out to be 
unstable, with the atomic planes adjacent to the cut al- 
ways slipping to either the perfect crystal or the intrinsic 
stacking fault configurations. Since the atomic density 
in the (111) plane is larger than in the (001) plane, the 
sp 2 hybridization of the atoms adjacent to the cut should 
be less favorable than the sp 3 , which would explain why 
platelets are not observed in the (111) plane. 

We now compare our model with the self-interstitial 
model. First, according to the results in Fig. 1, 
the activation energy per unit of area required to 
form the metastable stacking fault in our model is 
(3.42+0.75) J/m 2 = 0.26 eV/A 2 . This activation energy 
per unit of area should be compared with the activa- 
tion energy per unit of area required in the interstitial 
model. According to Goss et ai, the lowest energy to 
create a Frenkel pair at the platelet, for climb of the 
jog, and to diffuse the vacancy away from the platelet is 
about 8.8 eV, which would imply in an energy per unit of 
area of 1.39 eV/A 2 |25j- Therefore, the activation energy 
for the formation of the stacking fault in our model is 
about 5 times smaller than the corresponding activation 
energy in the model based on interstitials. Furthermore, 
we have performed preliminary calculations which indi- 
cate that the presence of nitrogen aggregates can fur- 



ther reduce the activation energy for the formation of 
the metastable stacking fault. The physical reason for 
this is the fact that since N is a trivalent impurity in 
diamond, less bonds have to be broken by the shearing 
process that creates the defect. Aside from a more favor- 
able energetics of formation, we believe that our model 
accounts for the optical properties of platelets in a more 
complete and straightforward way. Also, regarding sym- 
metry properties, the study by Goss et al. indicates that 
the interstitials can be arranged in a number of different 
structures, which have similar formation energy, not all 
of them exhibiting the asymmetry between the [110] and 
the [110] directions. In contrast, this asymmetry is es- 
sentially built-in in our model. In a nutshell, we propose 
a new model for platelets in diamond based on the shear- 
ing process of {001} planes that explains all experimental 
data available. 
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